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Investigations of molecular magnets are driven both by prospective applications in future storage technology 
or quantum computing as well as by fundamental questions. Nowadays numerical simulation techniques and 
computer capabilities make it possible to investigate spin Hamiltonians with realistic arrangements of local 
anisotropy tensors. In this contribution I will discuss the magnetic response of a small spin system with special 
emphasis on non-collinear alignments of the local anisotropy axes. 
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1. Introduction 

Since the early investigations of Mni2-acetate [ [1] single molecule magnets (SMM) are at the 
heart of the investigations of magnetic molecules worldwide. This is due to their properties which 
are governed by the anisotropy barrier, as there are slow relaxation of the magnetization as well 
as spin tunneling through the barrier [ [2j [3l IH [5] . From a chemical point of view it is striking that 
until recently no compound could be synthesized with a higher anisotropy barrier which moved 
hopes of easy application into the more distant future. Theoretical estimates pointed out that this 
might be due to very general reasons [ [HI E] . A recently synthesized manganese compound put an 
end to more than 20 years of search [ [5] . 

On the theory side numerically exact evaluations of spin Hamiltonians including anisotropic 
terms turned out to be limited to rather small systems such as for instance an antiferromagnetically 
coupled Ni4 compound [ (TUl [TTl [T^] , which was investigated in great detail or another but ferro- 
magnetically coupled Ni4 [[13|. With great numerical effort the exchange constants of Mni2-acetate 
could be determined [[Mj, but usually one resorts to models where the interacting spin system is 
replaced by one large spin in its effective ligand field (giant spin approximation). But thanks to 
the technical progress detailed numerical studies of larger molecules are affordable nowadays. First 
examples are given by the simulation of the new Mng compound with record anisotropy barrier [ 
[TS] and by the investigation of a MngFe compound [ilSj. Whereas Ref. [[15] restricts its parameter 
space to coUinear local anisotropy axes, Ref. [[16] considers non-collinear ones. 

In this contribution the recently developed procedures will be used to discuss the influence of 
non-collinear anisotropy axes on the magnetic response of a fictitious molecular compound. In order 
not to make the situation too complicated I will restrict the discussion to a triangular arrangement 
of equal spins of s = 1 with a molecular C3 symmetry. Then the magnetization will be discussed for 
three typical scenarios given by the possible ratios of exchange and anisotropy: | J| > | J| = 
and |J| < 
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2. Hamiltonian and evaluation 



2.1. Hamiltonian 

It turns out that many magnetic molecules can be well described by a Heisenberg model which 
takes the super exchange between the moments of the (typically 3d) ions into account. In this 
strong coupling limit anisotropic terms appear as perturbations and are thus often taken care of 
by perturbation theory. But for single molecule magnets, that might e.g. contain Mn(III) or V(III) 
ions, the situation might be reversed or anisotropy and exchange are at least of the same order. 
Therefore, I would like to treat both on the same footing right from the beginning. In the following 
Hamiltonian, 

ij i i 

the first term models the isotropic Heisenberg exchange interaction. A negative Jy corresponds to 
an antiferromagnetic coupling of spins at sites i and j. The second term models the local anisotropy 
tensors by their major principal axis which points along e.^. Depending on the sign of di this is an 
easy (di < 0) or hard {di > 0) axis. The neglected two other principal axes of the local d-tensor (e- 
terms) are usually much smaller. The last term provides the interaction with the applied magnetic 
field. Here gi is the local g-tensor, which might also be anisotropic. 




Figure 1. Configuration of the fictitious Cs-symmetric molecule discussed in this article. The 
remaining parameters of the model are the exchange coupling J, the strength of the local 
anisotropy d, and the angle of the local anisotropy axes with respect to the Cs axis. 

Since in this contribution the influence of the non-collinearity of the local anisotropy axes shall 
be discussed, the Hamiltonian will be restricted to a simplified spin system, that is depicted in 
Fig. [TJ Three spins s = 1 of a fictitious Ca-symmetric molecule will be considered that interact 
with one and the same exchange J. The local anisotropy axes are also related by the C3 symmetry, 
so that the anisotropy can be modeled by one strength d and one azimuthal angle z9. The y-tensors 
are assumed to be isotropic and of value 2. 

2.2. Evaluation 

Since the various parts of Hamiltonian ([1]) do not commute in general all eigenvalues and 
eigenvectors have to be calculated for each magnetic field B, i.e. for each strength and direction of 
the field. The magnetization is then evaluated using all eigenvectors. I would like to mention that it 
is also possible to obtain the magnetization without using eigenvectors by numerical differentiation 
of the energy eigenvalues Ei^{B), which are functions of B. 

Very often the investigated substance is only available as a powder. Then an orientational 
average has to be performed in order to be able to compare to experiments. Since this average 
cannot be performed analytically one sums over a finite set of special directions. In our evaluations 
we use for each absolute value B = \B\ special sets of directions given by points on the unit 
sphere. These sets are called Lebedev-Laikov grids [[IT]- The directions, which contribute with 
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various weights, and their total number are chosen such that the angular integration of polynomials 
x*^ X X z™, where k + l + m < 131, can be performed with a relative accuracy of 2~^^ [HZ]- In our 
evaluations we normally use a Lebedev-Laikov grid with 50 points. These points can be generated 
with publicly available software [[18]. Equivalently one may cover the unit sphere with the vertices 
of regular bodies such as an octahedron, a cube, or a dodecahedron. For practically all examples 
we evaluated in the past we can say that averaging over the directions given by the vertices of 
a dodecahedron is as accurate as averaging over a Lebedev-Laikov grid with 50 points. On the 
contrary, averaging only over a smaller number of directions like ice, ±y, and ±z is insufficient. 



3. Anisotropy versus coupling 

In the following main part the magnetic response of the spin system shown in Fig. [1] will 
be discussed for three scenarios: \J\ > \d\, \J\ = \d\, and |J| < \d\. In most cases the figures 
show powder averaged magnetizations for -d — 0, 10, 20, 30, 40, 50, 60, 70, 80, 90 degrees. Solid curves 
correspond to-d — 0, 10, 20, 30 degrees, dashed ones to -d ~ 40, 50, 60, 70 degrees, and dashed-dotted 
curves depict d = 80, 90 degrees. In order to pronounce the details all magnetization curves are 
shown at a low temperature of T = 0.1| J|. 

3.1. Strong coupling limit 




Figure 2. Magnetization for antiferromagnetic coupling J = —1 and easy axis anisotropy 
d — — 0.1|J| (l.h.s.) as well as hard axis anisotropy d = +0.1|J| (r.h.s.). ■& — 
0, 10, 20, 30, 40, 50, 60, 70, 80, 90 degrees. 




Figure 3. Magnetization for ferromagnetic coupling J = +1 and easy axis anisotropy d — — 0.1| J| 
(l.h.s.) as well as hard axis anisotropy d = +Q.1\J\ (r.h.s.). i9 = 0, 10, 20, 30, 40, 50, 60, 70, 80, 90 
degrees. 
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In the strong coupling limit the exchange energy is much bigger than the anisotropy energy, i.e. 
I J| > This is the case for many molecules, e.g. antiferromagnetically coupled iron or chromium 
rings [[111 [201 [111 112 [23 111 [25l IIS [23 [28]. These systems thus have in common that a description in 
terms of a plain Heisenberg Hamiltonian often provides a very good approximation [ [29l[30l[3Tll32j . 
If needed the anisotropic terms can then be included via perturbation theory. 

In Fig. [2] and [3] the magnetization curves are shown for combinations of anti-/ferromagnetic 
coupling J = =f1 and easy/hard axis anisotropy d = =f0.1| J|. In the antiferromagnetic case (Fig.Hj) 
the low-temperature magnetization curve is a staircase due to the successive ground-state level 
crossings in the growing magnetic field whereas in the ferromagnetic case the magnetization curve 
is practically given by the Brillouin function of the ground state spin S — 3. Since the anisotropy 
is small in magnitude its angular variation is not noticeable. 

3.2. Intermediate coupling 




qmbB/i]! gMBB/i:ii 



Figure 4. Magnetization for antiferromagnetic coupling J = —1 and easy axis anisotropy 
d = (l.h.s.) as well as hard axis anisotropy d = +1|J| (r.h.s.). i9 = 

0, 10, 20, 30, 40, 50, 60, 70, 80, 90 degrees. 
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Figure 5. Magnetization for ferromagnetic coupling J — +1 and easy axis anisotropy d = — 1| J| 
(l.h.s.) as well as hard axis anisotropy d = +1\J\ (r.h.s.). i9 = 0,10,20,30,40,50,60,70,80,90 
degrees. 

In the intermediate coupling regime exchange and anisotropy are of the same order. In the 
examples presented in Figs. [4] and [5] they have the same amplitude | J| = \d\. Looking at the figures 
one easily notices that the now larger anisotropy smears out the angle averaged magnetization 
curves. Nevertheless, the underlying structure - steps for the antiferromagnetic case and Brillouin 
function like for the ferromagnetic case - is still clearly visible. Interestingly the structure of the 
magnetization curve is most similar to the strong coupling limit for intermediate angles d around 
50°... 60°. 
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In detail the curves follow the following trend with varying At small 'd (solid curves), i.e. close 
to the collinear configuration, the influence of the anisotropy is strong since it acts cooperatively. 
Then for intermediate (dashed curves) the influence of the anisotropy is weakest compared to 
the strong coupling limit. This is understandable if one recalls that d = arccos(l/'\/3) = 54,73° 
is the angle of a perfect octahedral alignment of the anisotropy axes, i.e. the three axis have a 
pairwise angle of 90°. In this conflguration the effect of the anisotropic terms on multiplets of 
the Heisenberg Hamiltonian is either canceled or rather small. This explains why both at 50° and 
60° anisotropy effects are hardly visible. For ^ > 60° (dashed-dotted curves) the influence of the 
anisotropy rises again. 

3.3. Weak coupling limit 
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Figure 6. Magnetization for antiferromagnetic coupling J = —1 and easy axis anisotropy 
d — — 5|J| (l.h.s.) as well as hard axis anisotropy d = +5|J| (r.h.s.). i9 = 
0, 10, 20, 30, 40, 50, 60, 70, 80, 90 degrees. 




Figure 7. Magnetization for ferromagnetic coupling J — +1 and easy axis anisotropy d — — 5| J 
(l.h.s.) as well as hard axis anisotropy d = +5\J\ (r.h.s.). i9 = 0,10,20,30,40,50,60,70,80,90 
degrees. 

The weak coupling limit, for which the exchange energy is (much) smaller than the anisotropy 
energy, i.e. \J\ < \d\, is interesting since this case seems to be relevant for an improvement of single 
molecule magnets (SMM). The curves in Figs. [6] and [7] show that the effect of the anisotropy is 
- as expected - strong, and in the case of antiferromagnetic coupling (Fig. [6]) very strong. Here 
one can say, that a hard axis anisotropy together with antiferromagnetic coupling (r.h.s. of Fig. [6]) 
completely destroys any structure of the magnetization curve. A (meta-) stable magnetic ground 
state that is separated by an anisotropy barrier from its counterpart of opposite magnetization is 
either not created or the low-lying spectrum is so dense that the mechanism of a barrier breaks 
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down. Zero-field split multiplets do not exist any more and are thus an inadequate picture for such 
a situation. In the antifcrromagnetic case with easy axis (l.h.s. of Fig. [5]) one can at least for almost 
collinear alignments of the anisotropy axes obtain a non-vanishing ground state moment that is 
stabilized by anisotropy. 




T/I]l T/PI 

Figure 8. Magnetization as a function of temperature for g^B/\,J\ = 0.1 for ferromagnetic 
coupling J = -1-1 and easy axis anisotropy d — — 5|J| (l.h.s.) as well as hard axis anisotropy 
d = +5\J\ (r.h.s.). 1? = 0, 10, 20, 30, 40, 50, 60, 70, 80, 90 degrees. 



In the ferromagnetic case a hard-axis anisotropy (r.h.s. of Fig. [7]) also has the effect to weaken 
the magnetization compared to a Brillouin function. The interesting case is here again the case of 
easy-axis anisotropy (l.h.s. of Fig. [7]). If the magnetization for a small field, here gnB/\J\ = 0.1, 
is plotted against temperature (Fig. [8]) one again sees that a hard-axis anisotropy tends to smear 
out and weaken the magnetization (r.h.s. of Fig. [5]) whereas the curves appear sharper for easy- 
axis anisotropy (l.h.s. of Fig. [5]). It is worthwhile to investigate the energy spectrum of the latter 
configuration in some detail and compare it to the respective cases of intermediate and strong 
coupling. 



3.4. Possible SMM behavior 
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Figure 9. Energy eigenvalues at B = for J = 4-1 and d — — 0.1| J| (l.h.s.), d — — 1| J| (middle), 
and d = -5|J| (r.h.s.). i9 = 0. 



Regarding the question which of the cases is preferential for a good SMM the energy eigenvalues 
have to be studied. This is done for the zero-field case and J = +1 and d = — 0.1| J| (l.h.s. of Fig. [9]), 
d — — 1| J| (middle of Fig. [51), and d = — 5| J| (r.h.s. of Fig. (5]). The anisotropy axes are aligned in a 
collinear fashion, i.e. = 0, in order to act constructively. For a good SMM it would be necessary 
to obtain a low-lying zero-field split multiplet with both high total spin as well as high anisotropy 
barrier. In addition, these levels should be largely separated from higher lying levels. Looking at 
Fig. [5] one notices that in the strong coupling limit the multiplets can be well separated, but the 
ground-state multiplet is practically not split. In the opposite case of large easy- axis anisotropy 
the ground state multiplet is nicely split but energetically overlaps with other levels (but not 
too strong). Such a scenario was for instance discussed for the aforementioned Mng cluster [ I15j. 
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It turns out that the intermediate case, where exchange and easy-axis anisotropy are of similar 
magnitude provides the best chances for obtaining a good SMM: the ground-state multiplet is 
correctly zero-field split and the higher lying states are well separated. 

4. Summary 

In this article the influence of non-coUinear anisotropy axes on the magnetic response of a 
small spin cluster was discussed. The main result is that for good single molecule magnets ferro- 
magnetic coupling together with an anisotropy of coUinear easy axes that is of similar magnitude 
is preferential. Various other scenarios have been discussed. 
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